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■ Abstract 

ly-^ ', The radiative decays B*{D*) B{D)^ are investigated in the 

■ framework of light cone QCD sum rules. The transition amplitude and 
Qh! decay rates are estimated.lt is shown that our results on branching 



Oh! 



X 



ratios of D meson decays are in good agreement with the existing 
experimental data. 



1 Introduction 



One of the main goals of the future B meson and charm- r factories is a deeper 
and more comprehensive investigation of the B- and D-mcson physics. 

The radiative decays constitute an important for a comprehensive study 
the properties of the new meson states containing heavy quark. However, 
for interpretation of the experimental results we immediately deal with large 
distance effects. It is well known that QCD sum rules method [1] take into 
account such large distance effects and a powerfuU tool for investigating the 
properties of hadrons. Nowadays, in current literature, an alternative to 
" classical QCD sum rules method" , the QCD sum rules based on hght cone 
expansion is widely exploited. Main features of this version is that it is based 
on the approximate conformal non-perturbative invariance of QCD, and in- 
stead of many process-dependent non-perturbative parameters in classical 
"QCD sum rules", it involves a new universal non-perturbative parameter, 
namely the wave function [3]. This sum rules previously were succesful ap- 
plied for calculation the decay amplitude E ^ p7 [4] , the nucleon magnetic 
moment, Qt^nn and Qp^T^ couplings [5], form factors of semileptonic and ra- 
diative decays [6-9], the i^A'y* form factor [10], S — > p7 and D — > p7 decays 
[11,12], B*B7r and D*D7r coupling constants [13] etc. In this work we study 
the radiative B*{D*) — > B{D)^ decays in the framework of the hght cone 
QCD sum rules. Note that these decays were investigated Ref. [14,15], in 
the framework of classical QCD sum rules method. 

The paper is organized as follows. In section 2, we derive the sum rules 
which describes B*{D*) B{D)'y in the framework of the light-cone sum 
rules. In the last section we present the numerical analysis. 
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2 The Radiative B* B-f decay 



According to the general strategy of QCD sum rules, we want to obtain the 
transition amplitude for B* B'-f decay, by equating the representation of a 
suitable correlator function in hadronic and quark-gluon language. For this 
aim we consider the following correlator 

n^(p,9) = ^Jd^xe'P^<0\T[q{x)J^b{x),m^l5qm\0>F (1) 
in the external electromagnetic field 

Fa/six) = i{epqa - eaq^)e''^'' (2) 

Here q is the momentum and is the polarization vector of the electromag- 
netic field. The Lorentz decomposition of the correlator is 

= ieij,„ai3Pi^eaqpIl (3) 

Our main problem is to calculate 11 in Eq. (3). This problem can be solved in 
the Euclidean space where both, p'^ and p'^ = {p + qY are negative and large. 
In this deep Euclidean region, photon interact with the heavy quark per- 
turbatively. The various contributions to the correlator function Eq.(l) are 
depicted in Fiq. (1), where diagrams (a) and (b) represent the perturbative 
contributions, (c) is the quark condensate, (d) is the 5-dimensional operator, 
(e) is the photon interaction with soft quark hne (i.e. large distance effects), 
and (f) is the three particle high twist contributions. A part of calculations 
of these diagrams was performed in [12,14,15]. 

First, let us calculate the perturbative contributions, namely diagram (b). 
After standard calculation of the bare loop we have 

Hi = T^Nc / xdx / dy— — ■ — — — ^ 4 

47r^ Jo Jo m^x + m^x — p'^xxy — p'^xxy 



where A^c = 3 is the color factor, x = 1 — x ,y = 1 — y, p' = p + q, Qq and nia 
are the charge and the mass of the hght quarks. Note that the contribution 
of the diagram (a) can be obtained by making the following replacements 
in Eq. (4) : nib nT'a, &q ^ eg. The next step is to use the exponetial 
representation for the denominator 

1 1 fOO 

— = -. / ciaa"-^e-"^ 

{n-l)\Jo 

Then 

Hi = C xdx dylrubx + m„a;] H dae^"'-''+"'b^-P'''""''^^^ (5) 

4:71^ Jo Jo Jo 

Application of the double Borel operator ^(M^)^(M|) on IIi gives 

fti = ^ / dx^WribX + rriax] exp{- ^"'^ ^"^''^ (cti + 0-2)} (6) 

47r^ (7i + (72 io X XX 

where cri — ^ and a2 — In deriving Eq. (6) we use the definition 

BiM^-"" = 5{1 - aM'') (7) 
Now consider the combination 

-M-,^i)H-,^2)^ (8) 

St S S (7102 

which just gives the spectral density [16]. Using Eq. (6) and Eq. (8), for the 
spectral density, we get 



QaNc 



x-i 



pi{s,t) = -TT' i dx5{s - t)6{s - {mb + niaf )9{t - {mb + niaf ) 



xo 



47r2 
{nibX + ruax) 



X 

where the integration region is determined by the inequality 

sxx — {mix + mix) > (10) 



Using Eq. (9) for the spectral density, we get 

{ - (mt - ma)A(l, 7C, /) + mfeZn— — -} (11) 

1 -\- n — I — A(i, I) 

P2{s, t) = pi{a^ b, rua ^ nib, Qa ^ Qh) (12) 

2 2 

where k — I — ^ and 



A(l, K, I) = VI + k'^ + P-2k-21- 2kI (13) 
Finally for the perturbative part of the correlator we have 

n^^^ = ^ / ds- ^ -[(g„ - g,)(i - ^) + Q^in^] (u) 

Here we have neglected the mass of the light-quark. Applying Eq. (7), the 
Borel transformation, we get 

After simple calculation, for the double Borel transformed quark conden- 
sate contribution, we have: 

n^-^ = Q,< qq > J^e~^^ (16) 
and the 5-dimensional operator contribution is 



11 - Qb<qq> ^,^,e 



r mlml 1 1 . 1 ml 

For the calculation of the diagram corresponding to the propogation of 
the soft quark in the external electromagnetic field, we use the light cone 
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expansion for non-local operators. After contracting the b-quark line in Eq(l) 
we get 

U^ = i d X — — < 0\q{x)j^{mh+ Ii)j5q{x)\0 >f (18) 

J [Znyi — /e^ 

Using the identity T^^TaTs = 9nal5 + ic^pfsf^f^apf:) Eq. (18) can be written as 

d^x -f < Q\q{x)(jppq\Q >F (19) 

The leading twist-two contribution to this matrix element in the presence of 
the external electromagnetic field is defined as [4] : 

< q{x)api3q >f= Qq < qq> [ du(t){u)Faf}{ux) (20) 

J 

Here the function (t){u) is the photon wave function. The asymptotic form of 
this wave function is well known [4,17,18] to be, 

(t){u) = &iu{l - u) (21) 

where ^ is the magnetic susceptibility. 

The most general decomposition of the relevant matrix element, to the 
twist-4 accuracy, involves two new invariant functions (see for example [11,12]): 

< ^{x)apiiq >F = Qq<qq> {[ dux^(t)i{u)Fpfi{ux) 

J 

/ du4)2{u)[xi3Xr^Fpr,{ux) 
J 

- XpXr,Fpr,{ux) - X^Fpfi{ux)]} (22) 
In [11] it was shown that 

0i(«) = - u){2, - u) 

Mu) = -\{l-uf (23) 



So, for twist 2 and 4 contributions we get 

(f){u)du 



n*™+*^*^*^ = Qq<qq>{ 



ml — {p + uqY 



In order to perform the double Borel transformation we rewrite denominator 
in the following way: 

^6 ~ (p + ''J^qY = ~ (1 ~ ~ (p + ?)^'?^ 

and applying the Wick rotation 

~ (p + ^?)^ ^ ^& + (1 ~ ""jP^ + (p + ?)^'*^ 
Using the exponential representation for the denominator we get 

-^a,- 1 , 1 \ /M'2 1 



m|(M? + M|), 

The mass of B*{D*) and -B(£)) mesons are practically equal. So, it is natural 
to take Mf = M|, and introduce new Borel parameter such that Mf = 
M| = 2M2. In this case the theoretical part of the sum rules become 



n 



theor 



rr,.? ^2^2 _2 1 



b. m^mi ttIq 



+ Q,<m> - e-'""')lM^,t,{^) 

TfP' 1 

- 4(1 + ^)(0i(l/2)- 02(1/2)]^ (26) 



In the derivation of Eq. (26), we have subtracted the continuum and higher 
resonance states from the double spectral density. The details of this proce- 
dure are given in [13]. 

For constructing the sum rules we need the expression for the physical 
part as well. Saturating Eq.(l) by the lowest lying meson states, we have 



^phys ^ < mf.b\B* X B*\B^ X B\bt'y,q\0 > 
{m%. - p^){m\ - {p + qf) 



These matrix elements are defined as 



(27) 



< Q\q-i^h\B* > = e^fB*mB* (28) 

<B\U-f5q\0> = (29) 

nib 

<B*\Bj> = eappaPaqpe;e<^^h/mB (30) 

Here h is the dimensionless amplitude for the transition matrix element, 
, and rriB* are the polarization four- vector and the mass of the vector 
particle respectively, Jb is the leptonic decay constant and is the mass 
of the pseudoscalar particle, and e^- are the photon momentum and the 
polarization vector. Applying the double Borel transformation we get for the 
physical part of the sum rules 

h p-(m^*+m|)/2M2 

n^^^^ = fB^ruB^fBruB (31) 

Note that the contribution of three-particle twist-4 operators are very small 
[4], and thus wc neglect them (Fig. (1)). From Eqs. (26-30) we get the 
dimensionless coupling constant as 



h = 



^(ml.+ml)/2M^ 

fB*mB*fBmB 

{-^ / dse-'^[(Q, - Q,)(l - !^) + Qbln 
47r^ Jm? s 



-^r^ /-. rnhnl ml , 
+ <,,>e ^[Q,(l-^ + ^] 

+ (e-»?/"'-e— /"')[Q,0(i)M2 
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- 4Q,(l + ^)(0i(l/2)- 02(1/2))]} (32) 

3 Numerical Analysis of the Sum rules 

The main issue concerning of Eq. (32) is the determination of the dimen- 
sionless transition amphtude h. First, we give a summary of the parameters 
entering in the sum rules Eq. (32). The value of the magnetic susceptibility 
of the medium in the presence of external field was determined in [19,20] 



If we include the anomalous dimension of the current qcTa/sq, which is (—4/27) 
at II — m.}, scale, we get 

and 

< qq >= -{Q.26GeVf 

The leptonic decay constants /b(d) and fB*{D*) are known from two-point 
sum rules related to B{D) meson channels: = 0.14 (0.17) MeV [13,21] , 

fB*(D'*) = 0.16 (0.24) GeV [13,22,23,24], = 4.7 GeV , m„ = = 0, = 
(0.8 ± 0.2) GeV^, mB'(D') = 5. (2.007) GeV and mB(D) = 5. (1.864) GeV , 
and for continuum threshold we choose s^{sq) = 36 (6) GeV^. 
The value 4>i{u) — 02('w) is calculated in [11]: 
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From the asymptotic form of the photon wave function, given in Eq. (21), 
we get 



0(1/2) = 3/2x 



(33) 



Having fixed the input parameters, it is necessary to find a range of 
for which the sum rule is rehable.The lowest value of M^, according to QCD 
sum rules ideology, is determined by requirement that the power corrections 
are reasonably small. The upper bound is determined by the condition that 
the continuum and the higher states contributions remain under control. 

In Fig. 2 we presented the dependence of h on M^. From this figure it 
follows that the best stability region for /i is 6 GeV"^ < < 12 GeV"^ , and 
thus we obtain 



Note that the variation of the threshold value from 36 GeV^ to 40 GeV^ 
changes the result by few percents. We see that the sign of the amplitudes 
for and S+ are different. This is due to the fact that the main contribution 
to the theoretical part of the sum rules comes from the bare loop and the 
quark condensate in external field (last term in Eq. (32)). In S°(S+) cases, 
both contributions have negative (positive ) signs and therefore the sign of h 
is negative (positive). To get the dimensionless transition amplitude for the 
decay D* — > D7, it is sufficient to make the following replacements in Eq. 



-0.1 ±0.02 



fB+^fe+h = 0.2 ±0.02 



(34) 



(32): 



rUb rric, fs'-iB) fD*{D), Qb Q, 



(35) 



and Sob ^ sod 
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Performing the same calculations we get the best stability region for h as 
2 GeV^ < M2 < 4 GeV^ (Fig. (3)), and we find 



fDo^fuoh = 0.12 ±0.02 

fD+'fo+h = -0.04 ±0.01 (36) 

The signs of the transition amplitudes for Do and meson decays are 
different as in the B-meson case. 

Using the transiton amplitude "h", one can calculate the decay rates for 
B*{D*) B{D)^, which can be tested experimentally. For the decay width 
for radiative decay B*{D*) B{D)^, we get 

r(B*^So7) - 0.2SKeV (37) 
r(S+*^S+7) = 1.20KeV (38) 

and 

r(D*^Do7) = UAOKeV (39) 
r{D+*^D+-f) = 1.50 KeV (40) 

(41) 

In order to compare the theoretical results with experimental data for 
D-meson decays, we need the values of the D* — > Dn decays widths.We take 
these values from ref. [13]: 

r(L>*+ ^ D°7r+) = 32±b KeV (42) 
r(D*+^D+7r°) = 15±2 KeV (43) 
r(D*o ^ L)o7r°) = 22±2 KeV (44) 
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From eqs. (39-43), for the BR, we obtain 



39% 



BR{D 



1+ * 



3% 



(45) 



These results are in agreement with the CLEO data [25], which are 



We see that our predictions on branching ratio are in reasonable agree- 
ment with the experimental results. 
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BR{{Dl ^ L>o7) 



(36.4 ±2.3 ±3.3)% 
(1.1 ± 1.4 ± 1.6)% 
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Figure Captions 



Figure 1: Diagrams contributing to the correlation function 1. Solid lines 
represent quarks, wave lines external currents. 

Figure 2: The dependence of the transition amplitude h on the Borel pa- 
rameter square M^. Sohd hue corresponds and dashed hue to B'^ 
meson cases. 

Figure 3: The same as in Fig. 2 but for D meson case. 
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